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On the Improbability of a Random Distribution of the Stars in 

Space. 

By Kael Pearson, E.R.S., University College, London. 
(Received April 5, — Read April 21, 1910.) 

(1) The object of the present paper is to discuss certain statistical distri- 
butions of the stars, and to consider how far they are in accordance with our 
^experience so far as that extends. For direct tests we require far more 
numerous determinations of parallax than are at present available. 

The largest sample of parallaxes hitherto published is that due to Kapteyn,* 
but, for reasons elsewhere stated, I do not feel able to include this in my 
survey. Small series collected by Newcomb and determined by the Yale 
observers will be considered below. 

The main point of my investigation will, however, turn on certain parallax 
(constants, which frequently arise in stellar statistical investigations, and which 
are capable of determination without a knowledge of absolute parallax. 

I shall use the following notation: m = a stellar magnitude, D = an 
absolute stellar distance, it = cfD = a stellar parallax,f c being a constant, 
depending on the unit of measurement and the diameter of the earth's orbit. 
P = logio7r = logioc — logioD, is the logarithm of the parallax. A rule over 
any of these quantities will denote their mean values. cr m , ctd, cr , <j P will 
represent their standard deviations, or square roots of mean square deviations 
from the mean values. 

V will denote what is termed a coefficient of variation, or the ratio of a 
standard deviation to the corresponding mean value. 

Vd = cr D /D, Y„ = cTttItt, and cr P are the parallax constants which naturally 
arise in the course of stellar statistics, and which to some extent can be 
measured without the aid of absolute parallax determinations. 

The letter r will be used throughout for a correlation coefficient, the 
attached subscripts denoting the characteristic whose association is under 
discussion. 

(2) Let Li represent the light of a star ; then this light depends : (a) on the 
;size of the star ; (6) its intrinsic brightness ; (c) its distance from us ; and, 

* " The Parallaxes of 3650 Stars of different Galactic Latitudes/' J. C. Kapteyn and 
W. de Sitter, * Publications of the Astronomical Laboratory at Groningen/ No. 20, 1908. 

t I have retained tt for parallax in accordance with much astronomical usage, although 
it introduces some confusion when it = 3'14... occurs in the same equations, and still 
more if it appears in the form sin n. p and P are needed for many other purposes. In 
manuscript I invariably write n, but this never meets with the printer's approval. 
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possibly (cl), on the existence of interstellar absorption. Or, if X represent 
an absorption coefficient, we shall have 

or, logioLi == logioFi — 21ogi Di — \Di, where Fi depends solely on the 
constitution of the star. 

But if we measure light in terms of the stellace, i.e., the light, L , of a star 
of tenth-magnitude as unit, and use Pogson's value of the constant, we have 

logio Li = 4 — 0*4 mi -f logio Lo, 
and transferring to parallax, 

4 — 0*4 mi = logio Fi— logio L — 2 logi c + 2 logio tt—Xc/tt, 
or mi = J\ — 5 logio it -f 2*5 Xc/tt. 

Here f\ is entirely independent of the star's parallax ; it contains a constant,, 
10 + 2*5 logio Lo-f 5 logio c, which disappears when we consider deviations from 
mean values, and a quantity, 2*5 logio Fi, which depends solely on the intrinsic 
brightness and size of the star. With our present knowledge there is no- 
reason whatever for supposing any association between the distance of a star 
from the sun and its size or intrinsic brightness. We shall therefore suppose 
r n f = 0. Further, if there be an absorption term, we may suppose it so- 
small that for our first approximation we can put X = 0. Thus we write : 

?% = /i — 5 Pi, 
whence it follows that for any series of stars 

nil = fi — 5 Pi, 
o-m? — oy i 2 + 25cr Pi 3 „ 

It is accordingly clear that for any system of stars, the easily calculable 
a m must be greater than the value 5 a Fi obtained from any suggested 
theoretical distribution. We thus have a means of testing the possible non- 
applicability of various theoretical systems. 

Again, let there be a second magnitude, m 2 , arising from a light pheno- 
menon at the same or sensibly the same distance from us as m x . This can 
arise in the cases of variable stars and of double stars. Then we have 

^2= /a— 5 Pi, 

Cm 2 = <r/ a a + 25<r Pl 2 . 
Whence, if N be the number of stars dealt with, 

S (mi-mi) (m 2 ~m 2 ) S {fi-fi)\h~fr) , 9 . S (P-P) 2 
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Thus, if there be no association between the actual lights emitted, we 
should have 

Gr m 1 cr 'iruP^n l m 2 ==: ^O CTp , 

the left-hand side of which is determinable for any series by observation and 
calculation. This provides another method of ascertaining whether a 
theoretical distribution is applicable in the case of light phenomena at the 
same distance which are assumed to be independent. 
In general 



r AA = 



\Za m 2 — 25 <r P 2 \/<7 m2 2 — 25 o-p 2 ' 



which enables us to ascertain the actual relationship between light-giving 
capacities, if we can hypothesise a reasonable distribution for the stars under 
discussion. Thus, for example, we might, perhaps, assume to a first approxi- 
mation that double stars are scattered at random in space, and might inquire 
whether any relationship exists between the intrinsic brilliancies of the 
primary and secondary. 

(3) My proposed discussion will naturally break up into two parts. First, 
the consideration of what values these parallax constants take for theoretical 
stellar distributions ; and secondly, how far any system of stars observable 
by us, and for that very reason a selected system, is likely to be a random 
sample of such theoretical systems. 

I consider first some sample theoretical distributions. 

(4) Both Newcomb and Pickering have investigated in the first place the 
effect on parallax frequencies of the uniform distribution of stars through space. 
This appears an obviously reasonable suggestion, although, as we advance, 
weighty reasons will arise for modifying this standpoint. I take first : — 

The Values of Vd, V ff , and o-pfor a Uniform Distribution of the Stars 

throughout a Sphere of radius b. 

Let p be the density of the stars per unit volume, then 

•5 3 pD .= 4tirr 2 p rdr, or D = |6, 
Jo 

rb 

^7rb 3 po- D 2 = 4:7rr 2 p(r — lh) 2 dr, or <td 2 = -i^b 2 . 

Jo 

Hence : V D = Vii = °' 258 - 

Again 

§7rb d p7r = 4:7rr 2 p-dr, or tT = - j , 

* o 



$rt 



*b ( c 3 c \2 3 / c \2 

^irWpcTT 2 = I 4:7rr 2 p( z ~- - r ) dr 9 or a J = -l- 



Vpo-T 2 = 4:7rr 2 pl -j dr, or 
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Hence V, = V| = 0-577. 

Lastly, if ^ be the modulus = 0*43429, 



c 



fnrPpV = 4:7rr 2 p/jb log e ~-clr y 

Jo r 

c 4 itu 
and, integrating by parts, = A7rZ) 3 ^log e r+— q— A 6 * 

P = ,*(log.£ + |). 

Again P - P = ^ (log e i - log e t - ~ j 

= //,(&— log e r), 
where & = log e &-— |-. Thus we have 

4:7rr 2 pu 2 (Jc—log e r) 2 dr = 4:irfi 2 Q.og e r) 2 d(^r z ) — /n 2 ¥ 

o " Jo 



■§■ 7T& 3 /5 CTp 2 



47rJ' 



Hence, integrating twice by parts and evaluating the indeterminate terms, 

we find 

|-7r^a-p 2 = j'-n-^Kloge J) 2 -! log e 5 + f-F], 

or crp 2 = ^u, 2 , crp = -i^ = 0*1448. 

The special quantity generally needed, 25cr P 2 , is given by 

25 crp 2 = 2*778 fi 2 = 0*5240. 

Now the results reached are, supposing the stars dealt with to be uniformly 
distributed throughout space, of some interest. We see that they are quite 
independent of the radius b over which we extend the sphere of our 
observations. Further, we remark that <j m 2 > a m 2 , and a mim2 r mi7rh ought in 
all cases where we may reasonably assume a uniform distribution to come 
out greater than 0*524. 

(5) I will now proceed to modify these constants by supposing the stellar 
universe arranged in two spherical layers about the sun, each layer having 
uniform but individual density. Let Si be the density up to radius b h and S 2 
the density from b\ to b 2 . Then it follows that 

v - o-P5fifi i ^(h-h) 2 h s h d S 2 (S 1 ^S 2 ) -\i 
D I (V (Sx - S 2 ) + b 2 %) 2 J ' 

V - 0-^77 <f 1 , 4 (h — h) 2 hb 2 S 2 (Si — S 2 )~\ i 

and 25 .p 2 = 0*5240 + 25 ^ f 2 f* ^ ^^ . 

(bf($ 1 -~B 2 ) + b 2 d S 2 ) 2 

It is clear that to raise the values of V D , Y V9 and 25 cr P 2 , we must make S 2 
less than Si. In other words, in a general way these parallax constants will 
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be higher if the stellar universe is rarer at greater distances from us. But 
the changes are not very rapid. If we make b 2 = 2bi, then for Si = 2S 2 , we 
find 25<7p 2 only 0*748, and for Si = 9S 2 , a very improbable supposition, 
25crp 2 has only risen to 1*090. For b 2 = 4&i, the values for the same densities 
are 0*661 and 1*419. 

A dense nucleus round the sun appears the only way in which substantial 
change could be made in these parallax constants. 

(6) I do not stay to discuss these points, because another method of 
investigating the problem will throw light on other matters also. If we 
supposed our stellar universe to be spheroidal in form, or that we dealt with 
a spheroidal portion of it, we should clearly be approaching the same result, 
after a certain limit, as supposing the density to decrease with the distance 
from the sun. 

We need not reject such a system with the sun in the centre as necessarily 
idle, for we are only considering the nature of the changes which very rough 
approximations to the actual state of affairs will make in the parallax 
constants. 

In particular, the limit to the oblate spheroid may be looked upon as a 
flat stellar universe, and this " bunlike " universe is worth consideration. 
Again, the limit to a prolate spheroid is a columnar universe, and, the 
constants of such a columnar universe being determined, we can superpose 
two such* columnar universes, and obtain the parallax constants for a double 
drift system. 

The polar axis of our spheroid being 2b and the equatorial axis 2a, I find 
for the mean distance of the population from its centre 



D = |M i + 1 + ^ 



tan-^} =!*{-!, + OiV)!tan-a, 
J 1(1 4V) 2 V J 



V 
where 7] = (a 2 /b 2 —l)i, for the oblate spheroid; and 

1-e 2 /T + e 



D =#H 1 + 



8 - «, 

e v l — e 



lo gVr 



where e = (1— a 2 fb 2 f, for the prolate spheroid. Further 



for both. Hence Y D = ctd/D, can readily be found. 

In the case of the parallax we have for the oblate spheroid 

- — 3i tan" 1 ^ 
lr ~2b^~^T i 



Y7r= \/f^ lo ^ ot ^~\x 



X 

where cos % = bfa and rj = tan ^. 

E 2 
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For the prolate spheroid 



7p = -; — lO^g 
2 I € ° 



a/ 



1 — e' 



^V 



*~ V 3 



-v|rtani^ 



1, 



3{log e cot(i7r-J^)} 2 

where cos <fr = afh and e = sin i/r« 

These equations have been found by straightforward but somewhat 
laborious integration. The determinations of the mean, P, and standard 
deviation, <r Fi of the parallax logarithm are not so simple. 

We have 

c 



^ira 2 bpVj fju = 2 p 2 rrr cos 9 r d0 cl?* log e 

Jo Jo 



*tt r r 



2 




J 


o Jo 


yi 


2" 


47T 




Jo • 


4-7T 


rjrr 


3 

4 






r 



p log e - <rZ ( - r d ) cl sin 6, 







p ( r 3 loge — H - r 3 ) c£ sin 0. 



Let r = <zv, then ^ 2 = (l+^sin 2 ^)" 1 for the oblate spheroid, and 
sin# = r] 1 ^!/^ 2 — !. Change the integration to v, and we have 

on integrating by parts, where 



7 






and vanishes when nq = 0. Next, 

^TraPbpcr-p 2 = 2 



J 



p 2 Trr cos # r c£# cZr( /£ log e - — P ) , 

T / 



= ^irpfj? 



tjt 



cl (sin 9) cl (r 3 ) / log e -—v), 



where p = 4 -f 7. Whence 



$ira*bp<Tj? = ^Vl" I d(sin«)d(r3)(log c 8 ^-^ 



3 



J ^ (sin 0) 
Jo 

(fir 
cl (sin 6) 




6 

r 



= 4j£ ^ j*^ (sin ^ ^3 i oge 2 L + ^f log, I + f r 3 - p 2 r 3 ) . 



\ ~ r 3 
Put r = av again, and change to an integration with regard to v. We find 

2 = ^ f^og/ — + §log, ^^-^Wl^- 
if] j 1 \ av av J 
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Integrating by parts, the first term will be found to vanish between the 
limits, and we have 

or op 2 //* 2 = -1— 7 2 + I, 

where I is the integral. It will be found that 

1= _ 7 lo ge (l+^)-^[ 1 ^ /l - v2 lo ge ldv, 

b V ] b/a v v 

or 250-p 2 /^ 2 = 2-778-257 2 -257log e (l + 7j 2 ) 

JO^f'vjg! l fc 

To find the last integral we put 

n/i" 11 ^ 1 7 f 1 1, 1 , f 1 l-v/l^ 2 ! 1 7 

log e - dto = ~ log e - ^ — ^ log e - civ 

o/« J b/a >> b/a 



i(lo ge f)- r 



l-x/l-^ i 
logg - a-y . 



b/a V V 



b 

The latter integral is finite from to 1. The curve given by 

l-.-/l-i?a x 

« = log, - 

V V 

was plotted for the range to 1 and mechanically integrated. 
For the prolate spheroid by a similar process we have 



where y = 1 — L sin" 1 e , 

€ 

and *j?ltf = {-^-yiogia-^- ^-^ r^log^A;, 

e Ji -y 

or 

25^ 2 /yx 2 = 2-778-25 7 /2 -25 7 / log e (l-6 2 )- 50v/l ~ e2 r^EIiog^fe 
To find the integral, I take 

lOge V dV = loge VClv—\ log e V dfe 



-y 



I 



o ' Ji Ji v 

bL b .a\ f^V-v-c/--!.! 
-(log, 1 + 7 — ~ log^^y. 



^(io ge i-i +? )_r--zvi, 2 -i 

a\ a b] Ji v 






The subject of the latter integral always remains finite. The curve 

log e v was plotted and mechanically integrated up to the 



v— \z y v 2 — 1 



V 

values required. 
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The following table gives the chief stages in the determination of 25<r P 2 : 

Oblate Spheroid. 



Ratio of axes. 



*(io*£)' 



1 -— */l — v - i 1 j 
t log*, ~dv. 

bja v - v 



ajb 
a lb* 
ajb 
ajb 

ajb = 



1000 

100 

10 

n 

= 2 



6-6010 

■4-2986 
■2-0083 
•0 -8697 

•0-5207 



23 '85854 

10 -60380 

2 -65095 

'60347 

-24023 



'13572 
-13559 

-12877 
-09098 
-05937 



Prolate Spheroid. 



Ratio of axes. 



7 



— i LOge — — 1 + — J 

a\ a b /' 



&/«!,_ vV-1 



V 



loge vdv. 



hja 

b\a 
bja 
bja 
bja 



2 

3 

10 

100 

1000 



+ '3954 
+ -5648 

+ 0-8522 
+ -9844 
+ 0*9984 



-38629 

1 -29584 
14 -02585 

361 -51702 

5908 -75528 



-08902 
-16606 
-35055 
-48890 
-59240 



Substituting in the above equations, we obtain, after some arithmetic, the 
values of V Dj VU and 25 a v 2 . 
The following table results : — 





Batio of axes. 


7T. 

00 


* D* 


25<r p 2, 


' 


i 

! 

ajb = oo 
ajb = 1000 


0*391 


1-360* 


t « 


1 '771 


-391 


1 -360 | 


Oblate 
jpheroi 


o/5 = 100 
ajb = 10 


1 -378 
'914 


0*391 
0-385 


1 -354 I 
1 -210 


#/& = 3 


0*680 


0*343 


*842 ] 




ajb = 2 


0'622 


'309 


'680 




a[b = 1 


0*577 


'iOo 


*524 

i 




#/«■ ==• 2 , 


'628 


0*369 


0*753 


<D TJ 


&/« = 3 


0-703 


0*438 


1-104 




ft/a = 10 


1*085 


0*606 


2 *643 


O CD 

C .1 jimtwf 


i/a = 100 


2 *532 


*649 


4 *657 


PM Ph 


bja = 1000 


6-016 


0*650 


5 '148 




00 - 

i 


0-650 


5 *239f 

■ 



The limiting value of y is 1 -log, 2-- log, 1?, and again of 



<r P V = 2 -778-25 (l~log,2) 2 4 50 x '13572, 



The limiting value of e is 1 and of y' is 1, and again of 

<7 p 2 / M 2 = 2 *778 - 25 (1 + log, (1 - e- 2 ) ) - 50 ( log, 5-1^ 



27 *778. 
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I shall speak of 0*524 as the " normal " value of the constant of the stellar 
system, 25ap 2 , since it is the value it would take if the stars were distributed 
at random through space. It is clear that to suppose the distribution spheroidal 
and not spherical will allow of some increase from the normal value, but far 
ess for the " bunlike " than the columnar universe. It is desirable at this 
stage to compare these numbers with actual data, however slender. 

(7) Kapteyn's parallax determinations being excluded, I can only give 
results deduced from the 72 stars cited by Newcomb* as having the 
best determined parallaxes, and from 163 stars observed at Yale.f The zero 
and negative parallaxes must of necessity be excluded, because they lead to 
infinite distance D and infinite logarithm of parallax P. This reduced 
for Vd and 25a- P 2 the number of stars to 69 Newcomb and 122 Yale — 
admittedly most slender and imperfect material. 

We have, after some considerable arithmetic, 





v D * 


Y„ 


25<r P 2 . 


Yale stars 


0-258 
1-075 
1-017 


0-577 
0*924 
1*443 


*5240 

4*212 

3*111 





For Y„ I have, to see what effect it would produce, included the zero- 
parallax stars. The result is to separate more markedly the Newcomb and 
Yale series, and this was to be expected. For Newcomb excluded all but 
the best determined parallaxes, and so, many zero cases which would 
certainly be included by the Yale observers. For Vd the Newcomb and 
Yale stars give singularly accordant results, but wholly out of accord with 
the " normal " system. Nor is this to be wondered at. Parallax is subject 
to a very large error of observation, and our apparent an and av have not 
their true values, but values equal to v / ^d 2 + cte ] 2 and v / ov 2 + °e 2 2 > where 
<r El and oe 2 measure the variability in observational error in distance and 
parallax respectively. Now let us take a El = <td and a^ z = a n , i.e. assume 
the variability in error to be as great as the quantity to be determined. 
This cannot be true, because to admit it is to assert that parallax determina- 
tions are wholly unreliable, but it will serve to indicate that the maximum 
limit to error does not suffice to reduce our observations to the " normal " 
values. 

Dividing our apparent values by \/2 we find V D = 0*760 and 0719 

* "The Stars : a Study of the Universe,' 1902. 

t Chase, Smith, and Elkins, "Parallax Observations on 163 Stars," 'Trans. Yale Univ. 
Observ.,' vol. % pp. 1—207. 
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V„* = 0*653 and 1/020 for the Newcomb and Yale stars respectively. The 
reduction is considerable, but it will not bring the observations into accord 
with the extreme values of Yd for our spheroidal universe, still less with the 
theoretical " normal " values of Y D = 0*258, and Y v = 0*577. The theoretical 
and observed values for 25 a F 2 are still further out of agreement. 

It must be at once admitted that the argument against a uniform distribu- 
tion of the stars, based on these data, is not of much weight unless as 
confirming other results. In the first place we must suppose that we take 
a random sample of all the stars within a given distance (b) of the sun, 
i.e. with greater than a given parallax. In both cases the stars have been 
selected because they had big proper motions. "Newcomb's stars have been 
largely chosen because they had a parallax greater than 0*00, i.e. parallaxes 
only sensible in the third place of figures have been rejected, but we cannot 
assert that they are a random sample of all stars within a given sphere. Still 
the general agreement in Yd between the two independent series when we 
cut off all the zero (i.e. 0*00) parallaxes is remarkable, and carries some if not 
very great weight. 

Eeturning now to our fundamental equation 

a m 2 == cr/ + 25o- P 2 , 

we find for the 69 Newcomb stars a m 2 = 6*943, which is clearly greater than 
the maximum possible for 2 5 o-p 2 , i.e. 4*212, and leaves ay 2 a minimum value 
of 2*631. But if we take the Yale stars we find <r m 2 = 2*621, which is less 
than the apparent value, 3*111, of 25o- P 2 , so that if the latter were correct we 
should have a negative value for oy 2 . This can only be satisfactorily 
accounted for on the assumption that /and P are not independent in our 
selection. In other words, although there is probably no organic relationship 
between intrinsic brightness and distance, still whenever we take a series of 
stars by any form of selection, we introduce a spurious correlation, which 
will generally of course be negative, between intrinsic brightness and parallax ; 
that is to say, of the stars at a great distance — i.e. stars with small parallax — 
only those of great intrinsic brightness are admitted to our series.f The 

* For both series together (245 stars) Y^ = 1*361, or with maximum correction 0-962, 
quite incompatible with 0*577. 

t This question of a spurious correlation arising from the manner in which the 
observer selects parallax stars is not without considerable interest, and worth following 
up in the case of the Yale stars. Our general relation is m =/— 5P, which leads us at 
once to: oy 2 = o- m 2 + 10(r TO o- p r mp + 25o- p 2 and o-/(r p r fp = o- TO o- p r mp + 25 cr p 2 . Now we can 
directly find o- m , <r p , and r mT> . Omitting Algol we have for the 121 Yale stars with 
positive parallax: <r TO = 1*6190, o- p = 0*3540, r mv = -0*2179, whence we deduce oy = 2*1225 
and r fp = 0*6677. Now (see p. 48), we have : 

f~ a constant — 2*5 log 10 F 
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modification which this idea requires in our theory will be developed in the 
following sections. 

(8) In dealing with star series it is not unusual to include all stars up to 
a given magnitude, and counts have been made which are fairly reliable 
of all stars up to the 7*5 magnitude; indeed, we may be said to have 
practically complete lists of all the lucid stars ( < 6 mag.). It will be clear, 
then, that if we can modify our theory, so that we do not include all stars up 
to a given distance from the sun, but all stars up to a given magnitude, we 
shall be able to compare it much more satisfactorily with observation. Now, 
this process of cutting off from consideration all stars above a given magni- 
tude introduces a spurious correlation between magnitude and parallax, or 
between logarithm of parallax and the light emission. 

We reduce the " density " of the more distant layers of stars, because the 
fainter stars in these layers do not reach the minimum admissible magnitude. 

In order to reach any conclusion, it is, I think, needful to assume a law of 
distribution of the size and intrinsic brightness factor of the stars. 

Ee turning to our equation for the magnitude, we have (see p. 48), 
neglecting any absorption factor 

m — 10-f 2*5 logi L -{-51ogioC--2*5 logioF— 51ogi 7T 
= 10 + 2*5 logioLo— 2*5 Iogi F + 5//,log e D, 

where fi is again the modulus, and D is the distance. I shall put 
log e D — x, and instead of the/ of the previous work write briefly y for 

10 + 2-5 logioLo — 2-5 log 10 F. 

XXlUb 

m = y + hfix, 

and y will be a factor the variation of which depends solely on the size and 
intrinsic brightness of the individual star.* 

Now, we need some law of frequency to distribute the character y, and it 
is not unreasonable to assume it in the first place to follow the same 
distribution — the Gaussian — which has been adopted by many astronomers 
for other stellar characters. I shall consider this point again ; at present, it 
will be at least a guide to our theory, and indicate the possibility of 

where F measures the amount of light emitted at the star. Hence if / and P are 
positively correlated, F and P will be negatively correlated, or we find in the Yale data 
that the light emitted and the parallax of the star are substantially correlated. This, of 
course, is an indirect result of the further stars being only included in the list when they 
have relatively greater absolute brilliancy than the near stars. 

* x and y are used for the P and / of the earlier parts of this paper as more convenient 
for the following rather lengthy algebraic analysis. 
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extensions. We have now carefully to distinguish between the actual stars 
in the universe and those which will be found in our list, which we will 
suppose includes all, or at any rate a perfectly random sample of all, the 
stars up to magnitude m . . 

Let y be the mean value of y for all the actual stars, and suppose the 
frequency of y given by 

z = const, x er^-WI*** = const, x e"^ a »\ 

if we write 77 = y — yo for brevity. 

Assuming a uniform distribution of stars of density p, we have between 
DandD + SD 

4:7rpD 2 BD stars = 4?r i oD 3 81og 6 D = 4:irpe Zx dx. 

Hence, since the hypothesis is that the intrinsic brightness of the stars and 
their distance from the sun are independent, it follows that the number 
of stars having characters lying between y and-f-Sy, x and x-\-hx is of the 
form 

const, x g-fcW^V* dy dx = z x e~"^W +ar dy dx, say. 

[Now we are to take all possible values of x and y subject to the condition 
that y + bjjbx does not exceed m . Clearly our limit is to be given by 

1 
mo—yo = V + Sp®, or x = — (i/ — y), 

O fJb 

where vq is written for m —yo for brevity. Hence, if N" be the total number 
of the selected stars. 

r + oo rSZ^o-ti 



f + CO f5/A v u " 
e-hiW+te dx dy 
— 00 J — CO 



= ¥0 



+ «> 3 



e-hW+f^Cvo-y)^ 



J —00 



u4^=f)f + 'r»('^H 



_3_ 



^ —CO 



= -i^S/A* 10 n ) y/ 2ir <r y . 

This determines z = -== — <r ^ v° 10 ,* / and the constant of our 

\/2ir a y 

frequency distribution. We can now proceed to find the mean 7) and the 
standard deviation a^ of the emitted lights of the selected population of stars. 
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We have 



N"^ = z \ J 7]e~^ <T y z+Sx dx drj 



•oo -co 



3 a**' 



31 \ 5/i o/i / 



_0 O 3/ 3^A/ 3 8V 

= — e5/x\ 10 /a / —J- )A s /£ r K(Jy. 

6 \ DLL / 



Therefore, tj = — 

Similarly 



_3_ 
5/x 



y > 



and y = 3/0 + ?? = yo — 



3 



<T 



?, 



y 



V* 



C + CC f- (l/ Q — Y)) 



7 ? %-iWo- ! / 2 +3.r fa fly 



= i z o 



•CO —00 






•CO 



\ 5>L6 5yL6 / 



or 



S 2 + *7 2 = 



1 



9 



OY 



\Z2tt a-y <j y 2 + ^r-~ v^ 73 " ^ 



25/L6 2 



\/27r °V 

and accordingly cr^ = cr y . 

Thus, while the effect of cutting off the series at a given magnitude is to 
lower the mean light emitted, it makes no change whatever in the variability 
of light emitted. Neither the mean nor the standard deviation of the light 
emitted depends at all on the magnitude at which we close our series. 

I now turn to the mean and variability of the logarithm of the distance. 

1 



N# = z 



+ °° r--(^o—>?) 

OfJ. 



e W/<r y z % e Zx dx dy 



— oo -co 
r + co 



= i z o 



= i z o 



-co 



\Ofjb 61 



+ co 



e 



-*{*+&**)* I'* 



-co 



LV"°" i+ 2V 2 "V l,? + 



7- (Ty 



xez 



drj 



3 / + l<r^\ 



Hence 



1 13<j 2 1/ 3 cr 2 \ 1 1 _1 

)yU, 3 25yL6 J 5yL6\ O [JL j 6 (JL 6 



or 



1 1 

^ = — S (logio D) = ^ S (logio c-logio tt) = I (rn,Q-y) - 1 , 

P = tf + l ogl0 c— J-(mo-y). 
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This should be compared with the result on p. 50. It gives the modified 
form taken by the mean logarithm of the parallax. 

Further, since m = y + 5fix, 

m = y + h/jLX 

= 3/0 — ^ + Wo-?/o + - F -^ f 

O/jb D/Jb 3 

= mo -^ = m - 0*7238. 



o 



Or, we have the remarkable general proposition that : The mean oj any 
random selection of stars not including any over the m magnitude ought to be 
0*7238 of a magnitude below this limiting value, if the stars are distributed 
uniformly through space. 

I next pass to the standard deviation of x. We have 



i 



1ST (a x 2 + x 2 ) = z \ r^^VjV 2 dx dy 



•OO J — CO 



'K^rr'K^K^^-^ 



= 4-^o ebfx\ " ' 10 ix j | g 



— CO 



2 1 



drj. 



— q-TT— (^0 — ^) + | 
O 0/£ 

Whence, after some reduction, we find 

2 J. - J. 2 

But it is easy to see that o-i ogeD is identical with <r\ ogeirj or a-^j [jl of our 
previous notation. Thus 

25c7 P 2 = 2/^ + ^ 2 = 0-5240 -fa-/. 

We see accordingly that the quantity 25<r P 2 when we cut off at a given 
magnitude will for every magnitude be still constant, but it will exceed by 
the quantity <r y 2 the " normal " value on the assumption that a random 
sample has been taken of all stars. 

I now proceed to find the product moment of x and y with a view to 
determining their correlation. 



r+co 



N" (!Pxy + V x ) = z o e~i tfW+sa^ dx dy 



i 

57^0-1?) 



= i^o 



J — CO *> — CO 
+ 00 

e 

*J — CO 






Whence, after some reductions, and substituting for x and y, we find : 

1 



VxiJ " T^' 
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Eeverting to our original equation, 

m = y + 5fix, 
we have m = y + Ofjix; 

m—m = y—y + 5fi(x—x) 
.and <r m 2 = a^ + 10/jp*, + 25/xW 

_ 2 9 ~ 2 I ^ ^A 6 _J_ ^r 2 

C7y ^ C7y "1 1- C7y 

y 

9^ // 2 

= =^- = 0*5240. 
9 

Thus the effect of cutting off our sample at a given magnitude produces a 
most important result ; we replace the statement that a m 2 must always be 
greater than 0*5240 by a much more striking one : The standard deviation of 
the magnitudes of a series of stars cut off at a given magnitude is constant, 
■whatever that magnitude may be, and equal to the normal value of the 
25 a? 2 parallax constant. 

We have here a very simple test of the uniform distribution of the stars 
through space. 

I return to the result 

yxy — ,— ^r- "y • 

O Li 

"But r xy = f X y\<Ty<r x ~ — - Gy ~ —- 1 ~ — 



5^<r x 5<7 P v / l,+ 0-524/o- y 2 

There is thus always a negative correlation between x and y, which again means 
.■a positive correlation between the logarithm of the parallax and y, or, finally, a 
negative correlation between the logarithm of the parallax (or the parallax) 
and light actually emitted. In other words, cutting off at a given magnitude 
has introduced a spurious correlation between parallax and light emitted, 
between a star's distance from the sun and its size and intrinsic brightness. 
It is not probable that such a relationship actually exists in the totality of 
the stellar universe, but it should be apparently demonstrable in all actual 
collections of stars. For example, suppose the maximum light to be emitted 
by stars belonging to certain spectral (or colour) classes, then we should 
expect to find a correlation between parallax and spectral class of the above 
character. I think this " cutting off at a given magnitude " is most probably 
the source of the fairly high correlation (0*36) already noted by me as existing 
for the Yale observations between spectral class and parallax.* 

* < Monthly Notices, E.A.S./ vol. 68, p. 428, 1908. 
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Further, this cutting off at a given magnitude will also modify the natural 
relationship between magnitude and parallax. Taking 

m—m = y—y + 5/jl(x—x) 
and multiplying by (x — x), summing and dividing by the number of pairs in 
the usual way, we have 

jpmx — ^wJ^x^xm — "pxy ~r f^CTx 

— — — , /r * -I — «+- — rr * —~ . C 

bp* y + 9 + .V J ~ 9 ' 

Thus, p mx = 02413, and is a constant. 
Further, since <r m = f /*, we have 

11 1 1 



r 



xm 



3 " Vi+sk* V 1+ 5§£» 

Finally, the correlation of the logarithm of the parallax with the magnitude 

will be given by 

- 1 

a/1 ... a y 
V 0*5240 

From the value previously given for r xy = —r Fy we find : 

or we are able from observable quantities to determine the correlation of 
parallax with the actual light emitted.* It may be noted that when there is 
a high correlation between parallax and magnitude, a low correlation will be 
found between parallax and actual light emitted. 

(9) Before we proceed to compare the results of the preceding section 
with observation, it is fitting to consider how far they really involve the 
assumption made that the lights actually emitted by the stars obey a Gaussian 
frequency distribution. It is probably quite as legitimate to assume this as 
to suppose that stellar proper-motions follow the Gaussian law , but I wish to 
show that it is by no means needful to assume this in order to obtain the chief 
results reached in the previous section. I have taken the Gaussian curve in 
the first place because it has a certain sanction from astronomical usage, but 
our results are really much more general. It is open to anyone to argue that 
the universe is a dying system, that there are a majority of dark stars, or 
stars of small intrinsic brilliancy. If this be true, not only will the distri- 

* In the case of the Yale stars <r y 2 = oy 2 = 4*5051 and cr//0'5240 = 8*5975, and we 
find rvm = -0*3228 and n?/= 0*9464, against the observed values of -0*2179 and 0*6677. 
This is, of course, further evidence that the Yale series is not a random sample " up to a 
given magnitude " of stars distributed uniformly through space. 
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tmtion of intrinsic brilliancies give a skew frequency curve, but the distri- 
bution of the logarithms of the factors which measure actually emitted lights 
will do the same thing. Accordingly we need, to cover this case, a frequency 
distribution far more general than the Gaussian, and one that can represent 
all degrees of skewness, from a type like A to one like B. 





This is met by using the curve 



\P 



which ranges from r} =b to ^ = — oo, and which can take the form B if p be 
negative. 

I have discussed this curve at length in a memoir in the ' Philosophical 
Transactions/ vol. 186, A, p. 373. 

By supposing y or rj to have this type of distribution, we shall at any rate 
determine how far our results are peculiar to the form or symmetry of a 
Gaussian frequency. 

The form of the frequency surface will now be 



Z = 2o(l— 2f^ +a % 



and we have at once 



N= 36 



"& rr~( J 'o-->?} 



•00 J — QO 



\P 



1— ?[ ey^ +Sx dxd7j 



J —00 ^ ' 



Assume £ = (7- 



3 



have 



5 jul / 



( & -,) = x(l-|), so that X = »(,-#-), and we 



I. 



PI 



E" = -}# O 5M- 






Thus 



2 = 



6r(p+i) 



5 jU, 
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Next : 



N5 = 



^o 



•CO*' -co 



(v Q ~rj) i 



(l—^Yeyi+teTjdxdi) 



•g-zo 



X? 



+i 



o 



?e-ib[l-£)dS, 



v - r (jo 



= 5 1 



wv<H>? 







Now 7] is measured from the mean value, and therefore it is needful that* 
p + 1 = by. Hence 

by\ __3^ == _ 3 jp + 1 _1_ 

A/ 5/xX. 



^ = 6(1 



5yU, 7 2 



o 



OyC67 



Next 



N (or^ + T) 2 ) = 2 



~(^o-^) I r .\P 



.J 



— ,M ey^ +3x 'r] 2 dxdrj 



cc ^ — -co 






/» 00 



x^ 



+ 1 



\ X X 2 / 



Hence 



and 



2 • -2 _ £2 / ! 



9, 



<v + ^ 



x (?l 1)+ x^ 



X J rf \ OfJLy, 



tnati is 

Further, we may write 



'y 



°"i? ~ o" ijl 11 — 



5 fty; 
Ofxyj 



V 



Icr/Hl 3 ^ 



&fJL 



bfiy. 



Thus a^ and ?? differ from the values on the Gaussian hypothesis, i.e. 

3 • / 3 \ 

<x v and cr v 2 } by the introduction of the factor 1/ 1 — - — , which becomes 

J hfju » \ o/juyj 

unity when 7 ' = 00 . Still the two general propositions remain true, i.e. that 

* 'Phil. Trans.,' A, vol. 186, p. 374. 
t a-y 2 = {p + l)/y 2 ; see loc. cit., p. 373. 
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when we cut off our series at a given magnitude, the mean of the actual 
lights emitted is reduced by a constant quantity, and the variability of these 
lights is constant, whatever be that magnitude. 

Turning next to the mean and variability of the logarithm of the distance > 
we find 

~Nx = so l 1 """"/;' eyr)+P ' x ' xdxd V 

Hence # = il — _.— _+ (# + 1), 

0/6 O 5yL6 /&\ 

__ v 1 , 3 h 



5/a 3 25^ 2 __3 ' 

7 5/6 

yp 1 , 3 cr y 2 

C .. ■ o"*" OK ,,2 



5/4 3 25/4 2 - 



o/z/y 
which differs from the Gaussian frequency result solely by the presence of 

the factor 1 — - — in the denominator of the last term. 
Clearly we can write 

leading, as on p. 59, to the result 

V = i/A + logi c-4(^o-y). 

Now y is a constant whatever be m , for it equals yo—y. Hence the mean 
logarithm of the parallax, when we increase the magnitude at which we 
terminate our series, decreases by ^ (difference of cutting off magnitudes). 

Proceeding next to the value of <r x , we have 

•6 r^ (vo ~ y} \ 



~N(a x 2 + x 2 ) = Zq ^h) eyr)+ x2dxd V 



■00*' —00 



Substituting for tj in terms of £ integrating out and using the above value 
of x, we find after some rather lengthy reductions that 
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We have then, as on p. 60, 

25 <tp 2 = ~fi 2 + <r v 2 = 0'5240-f v* = 0-5240 + a/(l~~)~ 2 . 
9 \ o/jLy/ 

Thus, as before, 25crp 2 is a constant quite independent of m , but exceeding 
its " normal " value by a fixed amount. 

We may add here the value of the mean magnitude. As on p. 60 

m = f + 5ftx 

= m ~^ = m -0-7238. 
3 

Thus our general proposition of p. 60, that the mean magnitude of 
a system of stars up to a given magnitude is always less by 0*7238 than that 
magnitude, is established on the broader basis of non-Gaussian variability. 

Lastly, I take the product moment of x and v 

rb r— O0-7?) / ~i\p 



J —ccJ —00 



= lzo 5 " 



^ +j 8. (,„-(,) &2 r u. /, _ rv r 1 / _ h , &n_i 



\p 



+1 



/T-yt^- 6+ Tj-^^ 



whence, integrating on the same lines as before, by aid of the T-funetion, 
we find 

\ bay! 
B life crj = 0\, 2 + 1 /f^ + 2 5 yuV/ 

— rr 2 —9rr 2 4- J- /r 2 

= ^£- = 0-5240. 

Thus, precisely as before with the Gaussian law, we have obtained with this 
more general frequency the same result, namely, whatever be the magnitude 
at which we cut off, the standard deviation of the magnitudes of the series is 
constant and equal to the " normal " value of 25<7 P 2 for an unselected 
universe. 

We may conclude with the correlations 

L , 0'5240\-* 

rxm = V 1 + : 524oj = " rpn " 
rp/ = 1— r Pm 2 . 
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It follows therefore that all the results previously obtained for a Gaussian 
distribution hold for the much more general type of frequency now adopted. 
The only differences depend upon slight changes involving the constant factor 

— -\ 

Spy/ 

The distribution I have selected for the emitted lights enables us to 
suppose brilliant stars or even dark stars to be in a majority f It is sufficiently 
wide, I think, to cover any probable distribution of the intrinsic brilliancies 
of the stellar universe very closely. The main conclusions, however, that 

m = mo-0'7238, a m 2 = 0*5240, 

I suspect to be capable of simple proofs independent of any assumption as to 
frequency distribution. I have not stayed to seek for such proofs, but given 
those by which I have myself reached these results. They provide the mean 
magnitude and magnitude variability of any catalogue or random sample of 
stars up to a given magnitude, for example, of all the stars which could be 
seen with a telescope of given aperture — of course, on the assumption that 
the stars are distributed at random through space. They are thus funda- 
mental in determining whether such a distribution is the actual one or not. 
There are two points to be borne in mind here. First our discussion is 
perfectly independent of direction, we have only supposed that spherical 
shells about the sun as centre, if of equal volume, contain the same 
number of stars. These stars may or may not be clustered in individual 
shells towards one or more directions. Our propositions apply to the 
frequency at each distance, and have no relation to direction. 
Secondly, we have seen that 

<r m 2 = 0-5240 = 25 oy 2 , 

or takes the value that 25<rp 2 does for a uniform spherical stellar universe. 
I strongly suspect that this principle is general, namely, that the value of 
cr m 2 when we cut off at a given magnitude is equal to the value of the 
parallax constant, 25cr P 2 , for the total universe, whatever be its assumed form. 
If this be true, then the values given in the table on p. 54 will have 
additional meaning. They show the modifications possible in a m 2 when we 
suppose the stellar universe spheroidal rather than spherical in form — still 
maintaining a uniform distribution of stars. As I have indicated, an 
ellipsoidal rather than a spherical form for the stellar universe appears to 
raise the value of 25crp 2 and accords more closely, as we shall see, with 
experience, I hope later to discuss how far <r m 2 = 25cr P 2 is either approxi- 
xnately or absolutely true for other than spherical stellar distributions. 

(10) It remains now to consider how far the result? reached in the 

F 2 
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previous sections, (8) and (9), are consistent with observation. Luckily we 
possess very accurate counts of the stars up to magnitude 7*25. The best of 
these is probably that due to Prof. Pickering,* and I shalj. make use of that 
here. Even if in some of the categories the counts are not absolutely 
accurate, the divergencies will occur in the higher magnitudes, where they 
will, owing to the large numbers dealt with, make no impression on the 
values of the means or standard deviations. 

The table below gives the frequency up to each magnitude in the second 
column. In the third and fourth columns are given the actual and 
theoretical mean magnitudes up to the magnitude in the first column, the 
fifth and sixth columns give the square of the actual and theoretical 
standard deviations. The seventh and eighth columns give the differences of 
the observed and calculated values registered in the third to sixth columns. 
The means and standard deviations were calculated by the summation method 
discussed by Palin Elderton.j- 



1. 


2. 


3. 


4. 


5. 


6. 


7, 


8. 


Up to 
magnitude. 


Observed 
frequency. 


Mean. 


Sq. standard deviation. 


Differences. 


Observed. 


Theoretical. 


Observed. 


Theoretical. 


Am. 


Aow 2 . 


0*25 
0*75 
1-25 
1-75 
2-25 
2*75 
3*25 
3-75 
4*25 
4-75 
5-25 
5-75 
6*25 
6-75 
7'25 


5 

10 

18 

32 

58 

105 

193 

336 

589 

1067 

1972 

3562 

6284 

11004 

17955 


1 -440 
1*914 
2*409 

2 -874 

3 -357 
3-869 
4-388 
4-885 
5*368 
5-853 
6-297 


1-526 
2-026 
2-526 
3-026 
3*526 
4-026 
4-526 
5-026 
5-526 
6-026 
6-526 


0-423 
0-512 
-571 
0-619 
0-664 
0-689 
0-691 
0*688 
0-695 
0-711 
0-748 


-524 
0*524 
0*524 
0*524 
0-524 
-524 
0*524 
0-524 
0-524 
0-524 
0-524 


-0-086 
-0-112 
-0-117 
-0-152 
-0-169 
-0-157 
-0-188 
-0-141 
-0-158 
-0-173 
-0-229 


-o-ioi 

-0-012 
+ 0-047 
4-0*095 
+ 0-140 ! 
+ 0-165 i 
+ 0-167 i 
+ 0-164 
+ 0-171 
+ 0-187 
+ 0-224 



It will be obvious from this table that while our formulas <r m 2 = 0*524 and 
m = Wo— 0*724, express approximately the truth (never differing by a quarter of 
a magnitude from it), yet the observed values differ sensibly and almost regularly 
from the calculated values — the higher the limiting magnitude the greater is. 
the divergence. That is to say, the more stars we take into our consideration,. 

* "Distribution of the Stars," 'Annals of Harvard College Observatory,' vol. 48„ 
pp. 149—185. 

+ < Biometrika,' vol. 4, pp. 374—8. 
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the less accurate appears the assumption of the uniform distribution of the 
stars through space. I think we are fully justified on the basis of this table 
in stating that a random distribution of the stars through a spherical space is 
inconsistent with observed facts* Or, if such a random distribution exists, 
then there must be a correlation between distance and intrinsic brilliancy. 
Such an apparent correlation would be produced, of course, if light were to 
any extent absorbed in its transmission through space. The few stars of 
lower than the third magnitude may be reasonably supposed uniformly 
scattered through space, since they fairly closely satisfy the necessary con- 
ditions ; beyond this magnitude we have continuous deviations which cannot 
possibly be looked upon as the random or irregular variations from a chance 
scattering. 

(11) It may hardly be needful to bring further evidence that the stellar 
universe is not a random distribution through space of stars of varying size 
and intrinsic brilliancy. But it may still be of interest to illustrate one or 
two points from different material. I chose Burnham's ' Catalogue of 
Double Stars ' because the selection there is of a somewhat different character 
to a star count, and further because there are records of stars up to a much 
higher magnitude. Confining our attention to primaries only, we have the 



following table :- 



m 



Range of magnitude. 


Number of stars. 


w.* 


0* m- 


Up to 6 *2 

„ 8-2 
„ 10 -2 

j, JLo *o 


813 
4947 
9711 
9867 
9900 


6 '204 (6 -496) 
7*112 (7*496) 

7 -979 (9 -496) 
8*025 (11-496) 

8 -042 (14 -796) 


1-008 
1*087 
1*385 
1-517 
1 *596 



* Theoretical means are placed in brackets; of course, no stress must be laid on more than the 
first two or three entries, and the whole evidence is not comparable with that from total counts. 

It is clear that <r m 2 here starts a good deal higher than the star counts 
leave off. But the general tendency is the same, the mean does not increase 
fast enough, and the variability, instead of remaining constant, tends to rapidly 
increase. Very faint doubles may easily be overlooked, but, I think, this is 
far from explaining the whole evidence of the above table, which is in the 
same direction as we have seen in the total counts of lower magnitudes. 
Double stars even more markedly than the whole stellar system do not seem 
to obey the condition of uniform distribution through space. 

Thus from several directions we find evidence accumulating that the stars 
cannot, without regard to their size or intrinsic brilliancy, be considered as on 

* I owe this table to the kindness of Mr. P. F. Everitt. 
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the average uniformly distributed through space. The only way I can see to 
save such an hypothesis is either to assume a limited universe, not spherical 
in form, or to suppose light absorbed in transit. It will be interesting 
to inquire further into the effect of modifying the form of a uniformly dense 
stellar universe, and secondly to inquire whether when we confine our stars to 
one spectral class we get any difference in the general nature of the above 
deviations. 

In conclusion, I must acknowledge my very great indebtedness to 
Miss Julia Bell, M.A., for much aid in the laborious calculations, of which only 
the final results are tabled in this paper ; to Mr. E. Cunningham for a sugges- 
tion as to the best form in which to put the integral required in the case of 
the oblate spheroid universe ; and to Mr. C. H. Williams for mechanically 
integrating and tabling both oblate and prolate spheroid integrals. 



The Liberation of Helium from Minerals by the Action of Heat. 

By D. Orson Wood, B.Sc, A.E.C.Sc. 

(Communicated by Prof, the Hon. R. J. Strutt, F.R.S. Eeceived April 12, — Read 

April 28, 1910.) 

Introductory. 

It is well known that the helium contained in radio-active minerals can 
be liberated by the action of heat, but, hitherto, no exact quantitative 
experiments have been made on the subject, and it was thought desirable to 
conduct such experiments with a view :— 

(1) To ascertain whether heat may not be used as a substitute for chemical 
methods in the case of minerals and rocks for which those methods are not 
convenient,* and, if so, the temperatures necessary. 

(2) To throw, if possible, more light on the way in which the gas is retained 
within the mineral. 

A number of experiments have been made by Traversf on cleveite and 

* Chemical methods break down (a) when the helium-content is very small, so that 
large quantities have to be treated to liberate a measurable volume of gas ; (b) when no 
simple reaction is available, or when such a reaction gives rise to large volumes of other 
gases. Iron sulphide, which occurs in many geological horizons, is an example of this. 

+ Travers, 'Roy. Soc. Proc., 5 A, vol. 64, p. 140. The minerals were heated to. redness 
in a hard glass tube, and, therefore, to a temperature not exceeding 600° C. 



